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Abstract 



A two-dimensional body, exhibiting a slight rotational movement, moves in a rarefied 
medium of particles which collide with it in a perfectly elastic way. In previously re- 
alized investigations by the first two authors, Plakhov & Gouveia (2007, Nonlinearity, 
20), shapes of nonconvex bodies were sought which would maximize the braking force 
of the medium on their movement. Giving continuity to this study, new investigations 
have been undertaken which culminate in an outcome which represents a large qualitative 
advance relative to that which was achieved earlier. This result, now presented, consists 
of a two-dimensional shape which confers on the body a resistance which is very close to 
its theoretical supremum value. But its interest does not lie solely in the maximization of 
Newtonian resistance; on regarding its characteristics, other areas of application are seen 
to begin to appear which are thought to be capable of having great utility. The optimal 
shape which has been encountered resulted from numerical studies, thus it is the object 
of additional study of an analytical nature, where it proves some important properties 
which explain in great part its effectiveness. 

Keywords: body of maximal resistance, billiards, Newton's aerodynamic problem, 
retroreflector. 
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1 Introdution 

One area of investigation in contemporary mathematics is concerned with the search for 
shapes of bodies, within predefined classes, which permit the minimization or maximization 
of the resistance to which they are subjected when they move in rarefied media. The first 
problem of this nature goes back to the decade of the 1680s, a time when Isaac Newton 
studied a problem of minimum resistance for a specific class of convex bodies, which moved 
in media of infinitesimal particles, rarefied to such a degree that it was possible to discount 
any interaction between the particles, and in which the interaction of these with the bod- 
ies could be described as perfectly elastic collisions [I]. More recently we have witnessed 
important developments in this area with the broadening of study to new classes of bodies 
and to media with characteristics which are less restrictive: problems of resistance in non- 
symmetrical bodies [21 |3l IH |5l [6], in nonconvex bodies of single collisions [3 IH 13 [9] and 
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multiple collisions (lOl [HI US] , bodies of developable surfaces [9] , considering collisions with 
friction fT3] and in media with positive temperature [13] • However most studies which have 
been published have given special attention to classes of convex bodies. 

The convexity of a body is a sufficient condition for the resistance to be solely a function 
of singular collisions — all the particles collide at once with the body. This attribute allows 
us to considerably reduce the complexity of the problems dealt with. Even the various studies 
of classes of nonconvex bodies which have emerged, especially in the last decade, are based 
almost always on conditions that guarantee a single impact per particle — SI El [9] . Only 
very recently have there begun to emerge some studies supposing multiple reflections (see, 

e.g. [ini Ellin]). 

In the class of convex bodies, the problem is normally reduced to the minimization of 
Newton's functional — an analytical formula for the value of the resistance. But, in the 
context of nonconvex bodies, there is not any simple formula known for the calculation of the 
resistance. Even if it is extremely complex, in general, to deal analytically with problems of 
multiple collisions, for some specific problems of minimization the job has not been revealed 
to be particularly difficult, there even being some results already available [10] 111]. If, on the 
other hand, we consider the problem of maximization, then in this case the solution becomes 
trivial — for any dimension, it is enough that the front part of the body is orthogonal to the 
direction of the movement. 

And what if the body exhibits, besides its translational movement, a slight rotational 
movement? When we think of this kind of problem, we have in mind, for example, artificial 
satellites, of relatively low orbits, which do not possess any control system which could sta- 
bilize their orientation, or other devices in similar conditions. In this situation we imagine 
that, over its path, the device rotates slowly around itself. 

The problem of resistance minimization for rotating nonconvex two-dimensional bodies 
has already been studied in |12] [T5] : it was shown that the maximal reduction of resistance, 
as compared with the convex case, is approximately 1.22%. In its turn, the problem of max- 
imization of the average resistance of bodies in rotation is far from being trivial, in contrast 
with that which occurs when we deal with purely translational movement. This class of prob- 
lems was, therefore, the object of study of the work carried out by the authors in |16| [T7] : 
nonconvex shapes of bodies were investigated which would maximize the resistance that they 
would have to confront if they moved in rarefied media, and, simultaneously, exhibited a 
slight rotational movement. With the numerical study which was executed, various geomet- 
rical shapes were found which conferred on the bodies rather interesting values of resistance: 
but it was in later investigations, performed in the follow-up of this work, that the authors 
managed to arrive at the best of the results — a two-dimensional shape which confers on the 
body a resistance very near to its maximum theoretical limit. It is this latest result which 
now is presented here. 

The presentation of the work is organized in the following way. In section ]2] we begin by 
defining, for the two-dimensional case, the problem of maximization, which is the object of 
the present study. Then, in section ]3j we describe the numerical study which was realized in 
the tracking of the body of maximum resistance and we present the main original result of this 
study: a two-dimensional shape which maximizes Newtonian resistance. The two-dimensional 
shape is then the object of study in section ]4j where some important properties are shown 
which help to explain the value of resistance which it displays. In section ]5] we present 
the main conclusions of our study and include some notes on possible working directions to 
undertake in the future. Finally, in appendices [A] and [B] proofs of theorems [l] and ]2] are 



2 



provided. 



2 Definition of the problem for the two-dimensional case 

Consider a disc in slow and uniform rotation, moving in a direction parallel to its plane. We 
will designate the disc of radius r by Cr and its boundary by dCr- We then remove small 
pieces of the disc along its perimeter, in an e-neighborhood of dCr, with e £ of value 
arbitrarily small when compared with the value of r. We are thus left with a new body B 
defined by a subset of Cr and characterized by a certain roughness along all its perimeter. 
The essential question which we put is the following: up to what point can the resistance of a 
body B be augmented? More than getting to know the absolute value of this resistance, we 
are principally interested in learning what is the increase which can be obtained in relation to 
the smooth body (a perfectly circular contour, in this case), that is, learning the normalized 
value 

^ Resistance(.B) 
Resistance(Cr) 



It is possible, from the beginning, to know some important reference values for the normalized 
resistance: R{Cr) = 1 and the value of the resistance R{B) will have to be found between 
0.9878 ( |12l [15)) and 1.5. The value 1.5 will be hypothetically achieved if all the particles 
are reflected by the body with the velocity v+ (velocity with which the particles separate 
definitively from the body) opposite to the velocity of incidence v (velocity with which the 
particles strike the body for the first time), v"*" = — v, the situation in which the maximum 
momentum is transmitted to the body. It is also possible for us to know the resistance value 
of some elementary bodies of the type B. This is the case, for example, of discs with the 
contour entirely formed by rectangular indentations which are arbitrarily small or with the 
shape of rectangular isosceles triangles. As was proved in [TT], these bodies are associated 
with resistances, respectively, oi R = 1.25 and R = ^/2. 

Apart from being defined in the disc Cr, it is assumed that the body to be maximized is 
a connected set i? G R^, with piecewise smooth boundary dB. Therefore, let us consider a 
billiard in \ B. An infinitesimal particle moves freely, until, upon colliding with the body 
B, it suffers various reflections (one at least) at regular points of its boundary dB, ending up 
by resuming free movement which separates it definitively from the body. Denote by convS 
the convex hull of B. The particle intercepts the (9(convi?) contour twice: when it enters into 
the set convi? and in the moment that it leaves. L = |(9(convi?)| is considered the total length 
of the curve d{comrB), and the velocity of the particle is in the first and second moments 
of interception represented by v and v"*", and x and the respective points where they 
occur. As well, the angles which the vectors — v and v+ make with the outer normal vector 
to the section of 5(convi?) between the points x and x~^ are designated ip and ip'^ . They will 
be positive if they are defined in the anti-clockwise direction from the normal vector, and 
negative in the opposite case. With these directions, both (p as well as 99+ take values in the 
interval [— 7r/2, 7r/2]. 

Representing the cavities which characterize the contour of B by subsets fii, • • •) which 
in their total make up the set coiwB\B, the normalized resistance of the body B (equation M) 
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takes the following form (cf. [17J): 



R{B) 



\d{coiwB)\ 
\dCL\ 



L 



(2) 



being Lq = \d{convB)r]dB\ the length of the convex part of the contour dB, Li = |9(convi3)n 
with i = 1, 2, . . ., the size of the opening of the cavity Oj, and R{^i) the resistance of the 
normalized cavity rij, in relation to a smooth segment of unitary size, with 



1/2 fK/2 



1/2 J -it/2 



[l + COS (/?) — (/?)) COS tpdi^dx. 



(3) 



The function should be seen as the angle of departure of a particle which interacts with 
a cavity that has opening of unit size and is similar to Oj, with the similarity factor 
— see illustration of figure [T} 




Figure 1: Example of trajectory of a particle which interacts with a cavity ^ij 



From equation ([2]), we understand that the resistance of B can be seen as a weighted 
mean i^iLi/L = 1) of the resistances of the individual cavities which characterize all its 
boundary (including resistance of the convex part of the boundary), multiplied by a factor 
which relates the perimeters of the bodies convS and Cr- Thus, maximizing the resistance 
of the B body amounts to maximizing the perimeter of conv-B (|5(convi?)| < |9Cr|) and the 
individual resistances of the cavities rjj. 

Having found the optimal shape O*, which maximizes the functional (|3]), the body of 
maximum resistance B will be that whose boundary is formed only by the concatenation 
of small cavities with this shape. We can therefore restrict our problem to the sub-class of 
bodies B which have their boundary integrally covered with equal cavities, and in doing so 
admit, without any loss of generality, that each cavity fij occupies the place of a circle arc of 
size e <C r. As with Li = 2rsin(e/2r), the ratio between the perimeters takes the value 

|9(convS)| sin(e/2r) ^ . {^/f)'^ 

\dC^\ e/2r ^ 2^' ^' 

or that is, given a body i? of a boundary formed by cavities similar to J7, from ^ and Q, we 
conclude that the total resistance of the body will be equal to the resistance of the individual 
cavity fi, less a small fraction of this value, which can be neglected when e <C r, 

R{B) « R{n) - ^^Rm. (5) 
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Thus, our research has as its objective the finding of cavity shapes O which maximize the 
value of the functional O , whose limit we know to be found in the interval 



= 93, thus < I /_( /2 2 '^os ipdifdx = 1.5. 



3 Numerical study of the problem 

In the class of problems which we are studying, only for some shapes of 0, which are very 
elementary is it possible to derive an analytical formula of their resistance ([S]), as we saw in 
the rectangular and triangular shapes previously referred to. For somewhat more elaborate 
shapes, the analytical calculation becomes rapidly too complex, if not impossible, given the 
great difficulty in knowing the function ip~^ : [—1/2, 1/2] x [— 7r/2, 7r/2] — > [— 7r/2, 7r/2], which as 
we know, is intimately related to the format of the cavity 0,. Therefore, recourse to numerical 
computation emerges as the natural and inevitable approach in order to be able to investigate 
this class of problems. 

There have been developed various computational models which simulate the dynamics 
of billiard in the cavity. The algorithms of construction of these models, as well as the those 
responsible for the numerical calculation of the associated resistance, were implemented using 
the programming language C, given the computational effort involved (language C was created 
in 1972 by Dennis Ritchie; for its study we suggest, among the extensive documentation 
available, that which is the reference book of its language, written by Brian Kernighan and 
Dennis Ritchie himself, [IB])- The efficiency of the object code, generated by the compilers of 
C, allowed the numerical approximation of ^ to be made with a sufficiently elevated number 
of subdivisions of the intervals of integration — between some hundreds and various thousands 
(up to 5000). The results were, because of this, obtained with a precision which reached in 
some cases 10~^. This precision was controlled by observation of the difference between 
successive approximations of the resistance R which were obtained with the augmentation of 
the number of subdivisions. 

For the maximization for the resistance of the idealized models, there were used the global 
algorithms of optimization of the toolbox Genetic Algorithm and Direct Search" (version 
2.0.1 (R2006a), documented in [19]), a collection of functions which extends the optimization 
capacities of the MATLAB numerical computation system. The option for Genetic and Direct 
search methods is essentially owed to the fact that these do not require any information 
about the gradient of the objective function nor about derivatives of a higher order — as 
the analytical form of the resistance function is in general unknown (given that it depends 
on ip~^(x,ip)), this type of information, if it were necessary, would have to be obtained by 
numerical approximation, something which would greatly impede the optimization process. 
The MATLAB computation system (version 7.2 (R2006a)) was also chosen because it had 
functionalities which allowed it to be used for the objective function the subroutine compiled 
in C of resistance calculation, as well as the ip~^{x, (p) function invoked in itself. 



1 < supf^R{n) < 1.5 



(6) 
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3.1 "Double Parabola": a two-dimensional shape which maximizes resis- 
tance 

In the numerical study which the authors carried out in |16| [T7] . shapes of $7^ defined by 
continuous and piecewise differentiable / : [—1/2, 1/2] R"*" functions were sought for: 

flf = {{x,y) : -1/2 < X < 1/2, < y < f{x)} , (7) 

with the interval [—1/2, 1/2] x {0} being the opening. 

The search for the maximum resistance was begun in the class of continuous functions 
/ with derivative /' piecewise constant, broadening later to the study of classes of functions 
with the second derivative /" piecewise constant. In the first of the cases the contour of $7/ 
is a polygonal line, and in the second, a curve composed of parabolic arcs. Not having been 
able with these shapes to exceed the value of resistance R = 1.44772, we decided, in this 
new study, to extend the search to shapes different from those considered in ([7| . We studied 
shapes defined by functions x of y of the following form: 

= {(x, y):0<y<h, -g{y) <x< g{y)} , (8) 

where h > and g : [0, h] M.q is a continuous function with g{0) = 1/2 and g{h) = 0. 

The new problem of maximum resistance studied by us can therefore be formulated in the 
following way: 

To End supg R{Q^) in the continuous and piecewise differentiable functions 
g : [0, h] W^, such as g{0) = 1/2 and g{h) = 0, with h>0. 

Similarly to the study which was carried out for the sets J7y , in the search for shapes $7^, the 
functions g were considered piecewise linear and piecewise quadratic. If in the classes of linear 
functions it was not possible to achieve a gain in resistance relative to the results obtained for 
the sets Oj, in the quadratic functions the results exceeded the highest expectations: there 
was found a shape of cavity $7^ which presented the resistance R = 1.4965, a value very close 
to its theoretical limit of 1.5. There were also carried out some tests with polynomial functions 
of higher order or described by specific conical sections, but, not having verified any additional 
gain in the maximization of resistance, it was decided not to report the respective results. 
There therefore follows the description of the best result which was obtained, encountered in 
the class of quadratic functions x = ±g{y). 

The value of resistance of the sets J7^ were studied, just as defined in ([s]), in the class of 
quadratic functions 

ghAv) = + Py + 1/2, for < y < /i , 

where h > and a = (given that gh,f}{h) = 0). In the optimization of the curve, the 

two parameters of the configuration were made to vary: h, the height of the dVl^ curve, and 
/9, in its slope at the origin (^('(0)). In this class of functions the algorithms of optimization 
converge rapidly towards a very interesting result: the maximum resistance was reached with 
h = 1.4142 and j3 = 0.0000, and assumed the value R = 1.4965, that is, a value 49.65% above 
the resistance of the rectilinear segment. This result seems to us really interesting: 

i. it represents a considerable gain in the value of the resistance, relative to the best result 
obtained earlier (in jl6| I17j). which was situated 44.77% above the reference value; 
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ii. The corresponding set 17^ has a much more simple shape than that of set flf associated 
with the best earher result, since it is formed by two arcs of symmetrical parabolas, 
while the earlier one was made up of fourteen of these arcs; 

iii. this new resistance value is very near to its maximal theoretical limit, which, as is 
known, is found 50% above the value of reference; 

iv. The optimal parameters appear to assume value which give to the set a configuration 
with very special characteristics, as in what follows will be understood. 

Note that the optimal parameters appear to approximate the values h = \f2 = 1.41421 . . . 
and /? = 0. The following question can therefore be put: 

Are these not the exact values of the optimal parameters? 

The graphical representation of the function R{h, /?) through the level curves, figure [2] are 
effectively in concordance with this possibility — note that the level curves appear perfectly 
centered on the (\/2, 0) coordinates; marked on the figure by "+" . Note also, in figure [sj the 




(a) (b) 
Figure 2: Level curves of the R{h, (3) function. 



resistance graph R{h) for /? = 0, where it can equally be perceived that there is a surprising 
elevation of resistance when h \f2. Thus the resistance of the ^^^P cavity was numerically 
calculated with the exact values h = \f2 and /3 = 0, the result having confirmed the value 
1.49650. 

There is yet one more reason which suggests also an affirmative response to the formulated 
question. The shape of the set QP^^^ with h = \f2 and /3 = is a particular case with which is 
associated special characteristics which could justify the elevated value of resistance presented. 
The two sections of the shape are similar arcs of two parabolas with the common horizontal 
axis and concavities turned one towards the other — see figure |4j But the particularity of 
the configuration resides in the fact that the axis of the parabolas coincides with the line of 
entry of the cavity (axis of x), and that the focus of each one coincides with the vertex of the 
other. 

This shape of cavity appears to effectively deal with a very particular case. In contrast 
with what happened with all the other shapes which were studied, the integrand function 
of functional dsl) displays a rather smooth shape, presenting only a few small irregularities 
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R{h) 

1.5 1 ^ 1 




12 3 4 

h 



Figure 3: Resistance graphic R{h) for /? = 0. 




-0,5 -0,25 0,0 0,25 0,5 / 



(a) (b) 
Figure 4: (almost) Optimal 2D shape — the Double Parabola. 



for ip angles of little amplitude. Noting this characteristic, and taking into account that 
the integrand function almost does not depend on x, the resistance was calculated, for this 
shape in particular, using the rule of Simpson 1/3 in the integration in order to (/?. The 
double integration in the equation Q was thus numerically approximated by the following 
expression: 

R=-AxAip ^ Wk {l + cos {ip~^{xi,ipk) - ^Pk)) cos ifk, (9) 

j=Ar^/2-fl k=l 

with Wk = 2 for k odd and Wk = 1 for k even, Xi = — 1/2 + (i — l/2)Ax, Ax = l/N^, 
(fk = — vr/2 + kAif and Acp = n/Nip. Nx and N^p are the number of sub- intervals to consider 
in the integration of the variables x and ip (both even numbers), respectively, and Ax and Aip 
the increments for the correspondent discreet variables. Given that the shape J7^^>o presents 
horizontal symmetry, the first summation of the expression considers only the second half of 
the interval of integration of the variable x. 

In order to be easily referred to, this shape of cavity (figure |4^) will be, from here on, 
named simply ^^Double Parabola^\ Thus, in the context of this paper, the term "Double 
Parabola" should be always understood as the name of the cavity whose shape is described 
by two parabolas which, apart from being geometrically equal, find themselves "nested" in 
the particular position to which we have referred. 
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Since the resistance of the Double Parabola assumes a value which is very close to its 
theoretical limit, in a final attempt to achieve this limit, it was resolved to extend the study 
even further to other classes of functions g{y) which admit the Double Parabola as a particular 
case or which allow proximate configurations of this nearly optimal shape. In all these cases 
the best results were invariably obtained when the shape of the curves approximated the shape 
of the Double Parabola, without ever having overtaken the value R = 1.4965. It was begun 
by considering functions g{y) piecewise quadratic, including curves splines, without achieving 
interesting results; only for functions g{y) of 2 or 3 segments was it possible to approximate 
the resistance and the shape of the Double Parabola. Cubic and bi-quadratic functions g{y) 
were also considerecQ but in both cases the process of optimization brought them proximate 
to the curves of quadratic order, with the coefficients of greater order taking values which 
were almost zero. The problem was studied in the class of conical sections, considering, for 
lateral facets of the cavity, two symmetrical arcs either of an ellipse or of a hyperbole. Also 
in these cases the arcs assumed a shape very close to the arcs of the parabolas. 

The Double Parabola being the best shape encountered, and dealing with a nearly optimal 
shape, in the section which follows it is the object of deeper study, of an essentially analytical 
nature, where the reasons for its good performance are sought. 

4 Characterization of the reflections in the shape "Double 
Parabola" 

Each one of the illustrations of figure |5] shows, for the "Double Parabola", a concrete tra- 
jectory, obtained with our computational model. It is comforting to verify that, with the 
exception of one trajectory, in all the others the particle emerges from the cavity with a ve- 
locity which is nearly opposite to that which was its entry velocity. This is the "symptom" 
which unequivocally characterizes a cavity of optimal performance. Even in the case of the 
trajectory of the illustration (f), the direction of the exit velocity appears not to vary greatly 
from that of entry. 

If we analyze the five first illustrations, we may verify that there exists something in com- 
mon in the behavior of the particle: in describing the trajectory, the particle is always subject 
to three reflections. This appears to be a determinant characteristic for the approximation 
of the two angles of entry and exit. If, for example, we imagine three trajectories with prox- 
imate configurations, respectively, the trajectories (a), (b) and (c), but with the difference 
of not possessing the third reflection, the result would be completely different, as easily can 
be seen in the illustrations. Although this conviction is by nature essentially empirical, the 
results of the study which follow are heading in the direction of confirming that one very 
significant part of the "benign" trajectories — those in which the vectors velocity of entry 
and of exit are nearly parallel; we call them so because they represent positive contributions 
to the maximization of resistance — suppose exactly three reflections. 

We now will try to interpret another type of results obtained with our computational 
model, commencing with the graphical representation of the distribution of the pairs ((^, <p^) 
on the Cartesian plane — see figure [6| This graph was produced with 10.000 pairs of values 

^In the bi-quadratic curves, the point of interception of the trajectory of the particle with the boundary of 
the cavity is calculated by resolving an equation of the 4th degree. The roots of this equation were obtained 
numerically using the method described in [20] . The equations of inferior order were always resolved by utilising 
the known analytical formulas. 



9 



{x, if), generated by a random process of uniform distribution. 

The points concentrate themselves on the proximities of the diagonal = ip~^ , which 
revealing of good behavior on the part of the cavity. In addition, with these results it is 
shown that the response of the cavity deteriorates as (p approaches zero. Therefore, it begins 
to be understood that the "benign" trajectories have their origin essentially in entry angles 
of elevated amplitude. 

If we consider figure |6] there appears to exist an additional perturbation in the behavior 
of the cavity when the amplitude of the entry angle is inferior to about 20°, which means 
that some (99, (p~^) pairs become, in relation to the others, more dispersed and more distant 
from the diagonal (p~^ = ip. We have already called attention to the possible importance of 
the three reflections in the degree of approximation verified in the angles 99 and 99+. It occurs 




(a) X = 0.45, v? = 75°. (b) x = 0.45, 99 = 55°. (c) x = 0.45, ip = 35°. 




(d) x = 0.3, ip = 75°. (e) x = 0.0, ip = 35°. (f) x = 0.48, ip = 5°. 
Figure 5: Example of trajectories obtained with the computational model. 




Figure 6: Distribution of the y?"*") pairs on the Cartesian plane. 
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to us, therefore, to put the fohowing question: is it not precisely the number of reflections 
that, on differentiating themselves from the 3 occurrences, interfere so negatively with the 
behavior of the cavity? The investigations that follow will demonstrate, among other things, 
that our suspicion on this point has a basis. 

The following theorem says that for (p outside some interval {—ipQ^ipo), the number of 
reflections is always three. The proof is presented in appendix [A) 

Theorem 1. For ip entry angles superior (in absolute value) to ipo = arctan ( ^ ) — 19.47°, 



the number of reflections to which the particle is subjected in the interior of the Double 
Parabola cavity is always equal to three, and they occur alternately on the left and right 
faces of the cavity, no matter what the entry position may be. 

As a way to verify that the deductions which we have made are effectively in concordance 
with the numerical results of the computational model which was developed, we present one 
more graph, figure [7) produced with 10.000 pairs of (x, (p) values, generated randomly with 
uniform distribution. As can be observed in figure [7] all the trajectories with 4 or more 
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Figure 7: Distribution of the {ip, nr) pairs on the Cartesian plane, being nr the n° of reflec- 
tions. 



reflections, among the 10.000 considered, happened within the interval {—pQ,ipQ). Outside 
this interval (for \ip\ > ipo) the trajectories are always of three reflections. Additionally, we 
can verify that there isn't any trajectory with less than three reflections. This numerical 
evidence is confirmed by the following theorem: 

Theorem 2. Any particle which enters in the cavity Double Parabola describes a trajectory 
with a minimum of 3 reflections. 

The proof of theorem 2 is presented in appendix |Bj 

Of the conclusions which we arrived at we can immediately come to the following corollary: 
in trajectories with 4 or more reflections the angular difference \ip — 99"'' |, no matter how much 
bigger it may be, will never be superior to 2(/9o — 38.94°, a value which is much more inferior 
to the greatest angle which it is possible to form between two vectors (180°). The proof of this 
corollary is simple: as a trajectory of 4 or more reflections is always associated with a entry 
angle —990 < 95 < 95O) the exit angle will be situated necessarily in the same interval; taking 
into account the property of reversibility associated with the law of reflection which governs 
reflections, if just to be absurd we were to admit \p~^\ > po, on inverting the direction of the 
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particle, we would be in the position of having a trajectory of more than 3 reflections with a 
ip'^ entry angle situated outside the interval (— (^Qj'/'o)) which would enter into contradiction 
with the initial postulate. 
Summarizing: 

• There is verified a great predominance of trajectories with 3 reflections; 

• There are no trajectories of fewer than 3 reflections; 

• The critical angle ipo has the value lpq = arctan — 19.47°; 

• Outside the interval (— (/'Oj '■Po)-, ^ the trajectories are of 3 reflections; 

• In trajectories with 4 or more reflections, the angular difference is delimited by 2'Pq: 

5 Conclusion and future perspectives 

In the continuation of the study carried out previously by the authors in [HI [T7] , with the 
work now presented it has been possible to obtain an original result which appears to us to 
have great scope: the algorithms of optimization converged for a geometrical shape very close 
to the ideal shape — the Double Parabola. This concerns a form of roughness which confers a 
nearly maximal resistance (very close to the theoretical upper bound) to a disc which, not only 
travels in a translational movement but also rotates slowly around itself. In figure [8] one of 
these bodies is shown. Noting that the contour of the presented body is integrally formed by 




Figure 8: (almost) Optimal 2D body. 

42 cavities $7 with the shape of a Double Parabola, each one of which with a relative resistance 
of 1.49650, from ^ and ^ we conclude that R{B) = ^^^^^0^R{n) ^ 1.4951 is the total 
resistance of the body, a value 49.51% above the value of resistance of the corresponding disc 
of smooth contour (the smallest disc which includes the body). We know that if the body 
were formed by a sufficiently elevated number of these cavities, its resistance would even reach 
the value 1.4965, but the example presented is sufficient in order for us to understand how 
close we are to the known theoretical upper bound (50%). 

Although the value of resistance of the Double Parabola had been determined numerically, 
an analytical study was done in section [4j with the objective of consolidating the presented 
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results. We have managed to prove some important properties which help in the understand- 
ing of the elevated value of resistance which was obtained. We will try in the future to develop 
other theoretical studies which will allow us to consolidate this result even further. For ex- 
ample, an interesting open problem lies in delimiting the lag between the angles of entrance 
and exit for the trajectories of 3 reflections — for the others (trajectories with 4 or more 
reflections) we already know that \(p — 99"*"! < ~ 2 x 19.47°. 

The Double Parabola is effectively a result of great practical scope. Besides maximizing 
Newtonian resistance, it is exciting to verify that the potentialities of the Double Parabola 
shape found by us could also reveal themselves to be very interesting in other areas of practical 
interest. If we coat the interior part of the Double Parabola cavity with a polished "surface", 
the trajectory of the light in its interior will be described by the principles of geometrical 
optics, in particular rectilinear propagation of light, laws of reflection and reversibility of light. 
Thus, as the computational models which were developed by us to simulate the dynamic of 
billiards in the interior of each one of the shapes studied (where collisions of particles are 
considered perfectly elastic) are equally valid when the problem becomes of an optical nature, 
we can also look at 2D shape found by us in this new perspective. Given the characteristics of 
reflection which the Double Parabola shape presents we can rapidly conceive for it a natural 
propensity for being able to be used with success in the design of retroreflectors — see in |21] 
the exploratory study of its possible utilization in roadway signalization and the automobile 
industry. 

An incursion into the three-dimensional case, carried out in |2T], also showed that the 
Double Parabola is a shape of cavity which is very special. Our conviction of its effectiveness 
was strongly reinforced when we obtained the best result for the 3D case. This result was 
achieved with a cavity whose surface is the area swept by the movement of the Double Parabola 
curve in the direction perpendicular to its plane. The value of its resistance {R = 1.80) having 
been a little below the theoretical upper bound for the 3D case (i? = 2), to go beyond this 
value will be also an interesting challenge to consider in the future. 

For the 2D case we envision greater difficulty in going beyond the result which has already 
been reached — whether for the proximity which it has to the theoretical upper bound, or 
for the fact that we have already carried out, without success, a series of investigation with 
just this objective. 

A Proof of theorem [1] 

Consider a particle which enters into the cavity in (x,0), with the vector velocity forming 
an angle with the vertical axis, just as is found represented in the illustrations of figure |9] 
where we assume that the axis of symmetry of the cavity is the axis of the y and that its base 
^0^1 is placed on the axis x. In this way, the position of the particle at entry of the cavity 
assumes only values in the interval (—5, \) x {0}. 

Given the symmetry of the cavity in relation to its vertical axis, it will be enough to 
analyze its behavior for ipo < ip < 90°. The conclusions at which we arrive will be in this way 
equally valid for —90° < p < —po- 

We will analyze therefore in detail and separately each one of the sub-trajectories which 
compose all the trajectory described by the movement of the particle in the interior of the 
cavity. 

> 

Sub-trajectory BqBi 
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For (p > (fQ, we have the guarantee that the first reflection occurs in the parabolic curve of 

the left side of the cavity, just as can be easily deduced from the illustration (a). So that the 

particle collides with the left curve it is enough that the (p angle is superior to arctan(x/\/2), a 

magnitude which has as upper bound (fo = arctan(\/2/4). We thus have the initial trajectory 

of the cavity represented in illustration (a) by vector BqBi. 
> 

Sub-trajectory B1B2 

After colliding in Bi, in agreement with the law of reflection, the particle follows trajectory 

BiB-2. We prove that B1B2 has an ascendant path — illustration (a). We trace the straight 

line A1A2, segment, parallel to the initial trajectory of the particle BqBi, which passes through 

the focus of the left parabola (^i). Because of the focal property of this parabola, a particle 

> 

which takes the sub- trajectory A1A2, after reflection at A2, will follow a horizontal direction 

A2A3 (proceeding after its trajectory, after a new reflection, in the direction of the focus Aq 

of the second parabola). Upon the occurrence of the flrst reflection of the particle at Bi, a 

> 

point of the curve necessarily positioned below A2, the trajectory B1B2, which it will follow 
straight away, will be on an ascendant path, since the derivative g| of the curve at this point 
(Bi) is superior to the derivative in A2, where the trajectory followed was horizontal. 




Figure 9: Set of illustrations to the study of the trajectory of particles with entry angle 
</? > (^0 — 19.47°, in the cavity "Double Parabola". 
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Although we now know that B1B2 takes an ascendant path, nothing yet guarantees to us 
that the second reflection happens necessarily in the parabola of the right side. If we are able 
to verify that for ip = (pQ the second reflection is always on the right side, no what the entry 
position X is, therefore, logically, the same will happen for any value 99 > (/jq. This premise 



can be easily accepted with the help of illustration (a) of figure 10 for any value of > ipo, 
with the first reflection at a given point it is always possible to trace a trajectory for 
if = ^Pq which presents the first reflection at the same point Bi] the second reflection at the 
curve of the right side being for the case (p = ipQ, necessarily the same will happen for the 
trajectory with p > p^, since the angle of reflection will be less in this second case, just as is 
illustrated in the figure. Consequently, it will be enough for us to prove for ip = pQ, that the 
second reflection always occurs in the parabola of the right side, so that the same is proven 
for any which is the p > ipQ. 





(a) (b) 
Figure 10: Illustrations to the study of the second reflection. 



In illustration (b) of figure 10 the trajectory until the second reflection of a particle with 



the angle of entry p}Q [BqBi and B1B2) is shown. As one can conclude from the illustration, 
the B2 reflection only will happen on the curve of the left side if the a angle is less than a' . 
We have determined the value of the two angles. 

Being (xi, yi) the coordinates of the point Si, we will have tan(a') = — xi/(\/2 — yi), thus 

„' = arcan (^^^J^^) = -tan (^^^y^^) = -tan . (10) 

In order to arrive at the value of a we resolve the system of three equations, of unknown a, 
and /9, which are taken directly from the geometry of the actual figure 

a + (3 + 6 = -K 
f3 = Po + e 

arctan (jyi) + po + 9 = ^ 

The tangent line to the curve in Bi makes with the vertical an angle whose tangent has as its 
value the derivative g| of the curve at that point (in y = yi), where g| = 3^(3?/^ — 2) ~ 2?^- 
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Because of this, that angle emerges represented in the third of the equations by the magnitude 
arctan(yi/2). By resolving the system, the following result is obtained for a 

a = 930 + 2 arctan(yi/2) = arctan(\/2/4) + 2 arctan(yi/2). (11) 



Finally we prove that a > a', no matter what yi G (0, \/2) is. Of the equations (11) and 



(10), it will be equivalent to proving 

arctan(\/2/4) + 2 arctan(yi/2) > arctan((V2 + yi)/4). 

Given that < yi < \/2, both of the members of the inequality represent angles situated in 
the first quadrant of the trigonometrical circle. Because of this we can maintain the inequality 
for the tangent of the respective angles. Applying the tangent to both of the members, after 
effecting some trigonometrical simplifications, we arrive at the following relation 

1 4^2 - V2yj + 16yi ^ V2 + yi 
4 4-yf- V2yi 4 

Which, with additional algebraic simplifications, takes the form 



2/1 



(^V2yi + U + yl^/{4-yl-V2yi)>0. 



As < yi < \/2, it can easily be seen that both the numerator and the denominator of 
the fraction present in this latest inequality are positive magnitudes. Thus a > a', which 
contradicts the condition which was necessary so that the reflection B2 could occur on the 
curve of the left side, it therefore being proven, as we intended, that in no situation does the 



reflection B2 of illustration (b) of figure 10 happen on the curve of the left side. Logically, we 
can therefore conclude that the same happens for whatever is ip > ipo: the second reflection 

of the particle occurs always in the parabola of the right side. 

> 

Sub-trajectory B2B3 



We prove that the sub-trajectory has a descendant path — illustration (b) of figure|9 

Imagine, for this purpose, a sub-trajectory AqA2, parallel to B1B2 and which passes through 

> 

the focus Aq. The sub-trajectory A2A3 which will follow the reflection in ^2 — a point of the 

parabola on the right side situated below B2 — will be horizontal. The derivative of the curve 

in A2 being superior to the derivative value in B2, the sub-trajectory B2B3 will necessarily 

be of a descending nature. 

Even if we already know that the sub-trajectory is descendant, we have not yet shown 

that sub-trajectory in no situation conducts the particle directly to the exit of the cavity. 

Therefore follows the proof that the reflection B^ always occurs in a position superior to Aq 

— illustration (c) of figure [o] We trace A2B2, a segment of the horizontal straight line which 

passes through the point of reflection i?2- If the particle followed this trajectory, it would 

collide at the same point B2, but heads itself to Aq. Therefore, by the law of reflection, B3 

will have to be above Aq, since B1B2 makes an angle with the normal vector at the curve in 

B2 less than that formed by segment A2B2. 
> 

Sub-trajectory -B3-B4 

We will now show that the sub-trajectory which follows the reflection at B^ crosses the 
segment ^0^1 1 that is, directs itself to the outside of the cavity — illustration (d) of figure joj 
We trace, therefore, A2B3, a segment of horizontal straight line which passes through the 
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point of reflection B^. If the particle followed this trajectory, it would collide at and 

would head itself towards Ai. Therefore, by the law of reflection, the straight line where the 

> 

sub-trajectory B3B4 is placed will have necessarily to pass below Ai, since B2B3 makes an 

angle with the normal vector at the curve in B3 bigger than that formed by the segment 

A2B3. We have shown that the sub-trajectory crosses the axis of the x at a point situated to 

the left of Ai, but we have not yet shown that it occurs to the right of Aq. For that, we will 

have to prove that the third is the last of the reflections, that is, that in no situation does 

there occur a fourth reflection in the parabola of the left side. There follows this proof, of 

them all the most complex one. 

In order to prove that following the third reflection there occurs no other collision in the 

left parabola, we will show that a fourth collision — represented by B4 in the illustration 

(a) of figure 11 — has its origin always in an entry angle ip inferior to 920 • We will thus 





(a) (b) 
Figure 11: Illustrations to the study of a hypothetical fourth reflection. 



study the trajectory of the particle in the inverse order of its progression: we commence by 
admitting the existence of the sub-trajectory B^B^ of illustration (a) and we will analyze its 
implications in all the preceding trajectory. 



In illustration (a) of figure 11 are to be found represented the sub-trajectories B2B3 and 



B3B4. We begin by relating 02 with as, the angles which the vectors B2B3 and B-^B^, 
respectively, form with the vertical axis. For these purposes we resolved the system of three 
equations, of unknown 02, 03 and ^2, which are taken from the geometry of the figure]^ 

03 = P2 + as 
arctan (^ys) + (32 = ^ 
a2 + 03 + 132 = 

obtaining 

Q2 = 2arctan(y3/2) - as, 

in which arctan(y3/2) is the angle which the straight line tangent at the curve in B3 makes with 
the vertical — the inclination of the straight line tangent is given by g| \y=y^ = a|(iy^ ~ 
\) \y=yj. ~ \y^- turn, the angle as can be expressed in the following way: as = 



■^The variables denoted by Xi and yi, with i = 1, . . . , 4, represent the coordinates of i-th point of reflection, 
identified by Bi. 
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arctan((x3 - 2;4)/(y3 - 2/4)) = arctan((^y| - ^y|)/(y3 - 7/4)) = arctan((y3 + y4)/4), which 

permits us to write 02 in function only of the ordinates and ^4 of the extremes of the 

> 

vector B^Bi, 

02 = 2 arctan(j/3/2) - arctan((y3 + y4)/4). (12) 

In order to be able to prove what we intend — impossibility of occurrence of the reflection 
B4 — we need to find a lower bound for the ordinate of the position where each one of the 
four reflections occurs, or in other words, to determine {2/1, ?/2) 2/31 ylli just that 

yi > yl, 2/2 > 2/2, 2/3 > 2/3, 2/4 > yh V((p,x) g ((/po,^/2) x (-1/2,1/2). (13) 

It can easily be understood that y| = 0. We will therefore determine the other three lower 
bounds, commencing with 1/2 ■ 

We know that < ^4 < 7/3; therefore, from (112^ we take away that 



arctan(y3/2) < a2 < 2 arctan(y3/2) — arctan(y3/4). 

Being aware that 02 is situated in the first quadrant of the trigonometrical circle, we can 
maintain the inequalities for the tangent of the respective angles. After some algebraic sim- 
plifications, we obtain 

y3/2 < tan(a2) < 2/3 (l2 + 2/1) /16. (14) 

The equation of the straight line which connects B2 to B3 takes the form x = m{y — 2/3) + X3, 
with m = tan(a2) and X3 = jy| — As we are interested in finding the ordinate of the point 
of interception of this straight line with the parabolic curve situated on the right side, with 
equation x = —\y'^ + we have to resolve the equation of second degree, in the variable 
y, which results in the elimination of the variable x by combination of the two previous 
equations. The ordinate y2, of the second reflection, being the positive root of the equation, 
takes the form 2/2 = —2m + ^^ Am? + ^my^ — y| + 4 . 

The magnitude y2 is expressed in function of two variables, m and 2/3, which as we know 
assume only positive values. So as to accept more easily the deductions which we are going 
to make in the tracking of 7/21 will imagine, without any loss of generality, that yz is a 
fixed value. We begin by showing that the derivative of 7/2 in order to the variable m. 



2 {2m + 2/3 - v^4m2+4m?/3-2/i + 4 
dm y^4r7?T4my55^^y|T4 



(15) 



has a negative value for whatever value of y^ is. As 2/3 < \/2, inevitably y| < 4, thus the two 
radicands (4m^ + 4:myz — y| + 4) present in the equation ( |15| ) have always a positive value. 
The restriction 7/3 < \/2 allows us still to successively deduce the following inequalities 

yi < 2 4^ 2y| < 4 <^ y| < 4 - y| <^ 4m'^ + Amy^ + yj < An? + 4my3 + 4 - y| 



<^ (2m + y3) < 4m + 4my3 - y| + 4 <^4> 2m + y3 < y Aw? + 4my3 - y| + 4 . 
This last inequality confirms that 4^ < 0, for whatever ys may be. In this way, the value y2 is 



so much less the greater is the value of m. As is examined in ( 14 ), m < M = y3 (12 + y|) /16, 
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thus y2 > -2M + y^4:M'^ + AMy^ - y| + 4 . Substituting M, there is obtained, after some 
simphfications, 

3 1,1 



1/2 > /(ys), with fivs) = "2^3 - gyf + ^^J272yi + 40y| + + 256 . (16) 
In order to find the minimum value of /(ys) we begin by computing its derivative: 

_6_ ^ 272^3 + SOyj + 3yl - (12 + 3y|)^/2^2yi + 40yl + j/g + 256 
dys 8\/272y|+40y| + y« + 256 

The radicands being clearly positive, we only have to concern ourselves with the numerator 
of the fraction. To find the roots of the derivative function 3^/(2/3) is equivalent because of 
this to resolving the equation 

(272^3 + 80yi + Syf)' = (l2 + 3ylf (272yi + 40yl + yl + 256) , 
which can be simplified in the following: 

2304 - 1024y| - 992y| - 160y^ - 3yl = 0. 

This polynomial equation has only one real positive root, of the value 2/3 = |\/— 51 + 6\/79 ~ 
1.017, signifying that /{ys) has a global minimum in 7/3, because, as we show in what follows, 
^jpfius) > and the function does not presents other points of stationarity. 

We thus show that j-^fius) > 0, with 



d2 ^ _ 18240y| + 2960yf + IBOyf + 3y^° + 34816 + 30720y| 
dy| 4(272y| + 40y| + + 256) § 

(816yi + 120y| + 3y| + 768y3) ^272^1 + 40y| + yf + 256 



(17) 



4(272y| + 40y| + yf + 256); 
We show that ^^/(ys) > is equivalent to showing that the numerator of the first fraction is 



superior to the numerator of the second, in (17). After elevating the two terms to the square, 
we arrive at the inequality 

-3yl^ + 36y;^^ + 9464^;^^ + 191616y^° + 1514112^^ + 5817344yf 
+13535232y^ + 15532032y| + 9469952 > 0. 

We easily prove the veracity of this relation, given that we have a unique negative term 
(—3^3^) which, for example, is inferior in absolute value to the constant term (9469952), 
ys < V2 3^3^ < 768 < 9469952. The proof is thus complete that /(ys) has a global 



minimum in ys, of the value /(y3) = |\/— 51 + 6\/79- 



Accordingly, from (16), we finally conclude that 



y2 >y2 = nV + 6\/79~ 1.356. (18) 
9 

A lower bound is then found for the height of the second reflection B2 (illustration (a) of 



figure 11). We now determine ^3 , a lower bound for the height of the third reflection 
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So that the reflection B2 occurs in the parabola of the right side it is necessary that the 
angle 02 is greater than the angle formed between the vertical axis and the segment of straight 
line which unites -B3 with the superior vertex of the cavity, 

a2 > arctan (^^) = arctan (4^) = -ctan (^^) . (19) 



This inequality, in conjunction with the second relation of inequality of (14), allows us to 
write 



/2 + y3)/4 < tan(a2) < ys (l2 + yl)/l6 ^ {V2 + ys)/^ < ys (12 + yl)/16, 
from which results the inequality 

yl + 8y3 + 4V2 > 0. (20) 



As the polynomial of the left hand-side of (20) has a positive derivative and admits a 
unique real root, we immediately conclude that it constitutes an inferior limit for 7/3, this 
limit being 

-.1 ^ 
y* = - (54^/2 + dVMd) ' - 8 (54V2 + 6\/546) ' ~ 0.670. (21) 



3 

It is left to us to determine yl, a lower bound for the value of yi — the ordinate where the 



first reflection occurs. To this end we resort to illustration (b) of figure 11 which gives us a 
more detailed representation of the part of the cavity where the first two reflections occur, 
Bi and i?2- The scheme presented was constructed counting that the first reflection (Bi) 
occurs at a point which is more elevated than that of the third reflection (-B3). This is in 
fact the situation. This itself can be proven by showing that 02 is always smaller than the 
angle formed between the normal vector at the curve in B2 and the vertical axis, or that is 
02 < f — arctan (I2/2) • 

Taking, in ( |14[ ), at the upper limit of tan(a2) and having in mind that y < \/2, we build 
the following sequence of inequalities which proves what is intended: 

~51.06° ~54.74° 



/2/3(l2 + yi)\ ^ (7V2\ vr f V2\ tt fl 

a2 < arctan I — j ^ arctan I I 2 ~ ''^ \ ~2 / ^ 2^ ~ arctan I -y2 

We now try to find y\ . We can define y\ as being the ordinate of the point of interception 

of the left parabola with the semi-straight line with its origin at point B2 , positioning as low 

as possible (7/2 = 2/2)' with equal slope to the largest value permitted for the slope of the 

> 

trajectory which preceded B2 {B1B2). The equation of the straight line which connects Bi 
to B2 takes the form x = m{y — 2/2) + ^2, with m = tan(ai) and X2 = —\y2 + 5- As we are 
interested in finding the point of intersection of this straight line with the parabolic curve 
situated on the left side, with equation x = — ^ , we will have to resolve the equation 
of the second degree, in the variable y, which results in the elimination of the variable x by 
combination of the two previous equations. Although we have two positive real roots, we are 
only interested in the smaller of the two, which takes the form 



yi = 2m- J Am? - 4my2 - 2/1 + ^ • (22) 
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As we said, if we do y2 = 1/2 and place the maximum slope to the straight line, which in 
the previous equations is equivalent to considering m minimum, we obtain yi = yl. Being 
m = tan(ai), we should determine the value of qi through the system of equations 

ai = 202 + 02 
arctan (52/2) + 02 + 6*2 = f 



which is taken from the geometry of illustration (b) of figure 11 Is obtained 



ai = vr — 2 arctan(y2/2) — 02- (23) 



From this latest equality, from ( |14[ ), and given that y2 < \/2 and y^ < \/2, we deduce that 
ai > vr — 2 arctan(y2/2) — arctan(?/3 (l2 + y|) /16) > tt — 2 arctan(l/\/2) — arctan(7\/2/8), thus 



m = tan(Qi) > tan(7r - 2 arctan(l/\/2) - arctan(7\/2/8)) = |§\/2. If in ([22]) we do m = §\/2 



and 2/2 = ?/2 then obtain the inferior limit for yi 



yl = ^\/2 - ^ Y 444498 -33120\/2y -51 + 6V79 -38400V79 ~ 1.274. (24) 
Summarizing, (y^, y^, y|) ~ (1.274,1.356,0.670,0). 



With the help of illustration (b) of figure 11 we will, finally, analyze the entry angle (f of 
the particle. With the system of equations (^yi = g| \y=yi) 

201 + 93 + ai = vr 
arctan {\yi) + + 6*1 = | 



and with equalities (23) and (12) we obtain, successively, 
99 = ai — 2 arctan(yi/2), 

if = vr — 2 arctan(yi/2) — 2 arctan(y2/2) — 02, (25) 
(/3 = vr — 2 arctan(yi/2) — 2 arctan(y2/2) — 2 arctan(y3/2) + arctan((y3 + y4)/4). 



Taking (25), from (|19|) we deduce that 99 < vr — 2 arctan (yi/2) — 2 arctan (7/2/2) 



arctan((v2 + y3)/4). In agreement with the definitions (13) and with the values found in 



(18), (21) and (24), we can conclude that 

99 < vr - 2arctan(y;f/2) - 2 arctan(y^/2) - arctan((\/2 + y|)/4) ~ 19.18°, 

or that is 

93 < 990 ^ 19.47°. 

With this we can finally conclude that it is impossible to have a fourth reflection, since for 
this to happen the particle would have to have entered in the cavity with an angle 93 < 990, 
as we have just finished showing — something which would contradict our initial imposition, 
99 > 999. As the cavity presents symmetry in relation to its central vertical axis, the conclusion 
to which we have arrived is equally valid for 99 < — 93O) thus being proven that which we 
intended (theorem [T]) : 

For 1 99 1 > 990, there always occur three reflections, alternatively on the left and 
right facets of the Double Parabola cavity. 
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B Proof of theorem [2] 



For \lp\ > the statement of theorem [2] is already proved. It remains to consider the case 
—LpQ<(p<ipQ, but given the symmetry of the cavity we will only need to study the interval 
< < (/?o. 

We will base our proof on some of the deductions which we made in appendix |A] the 
illustrations of figure [9] being especially useful to us. We will also assume to be truthful the 
following premise: "If the second reflection happens on the same facet of the cavity where 
the first reflection occurred, there will necessarily be a third reflection" . We exempt ourselves 
from proving this principle because it appears evident to us. 

For Q < ip < (fQ the first reflection can just as well occur on the left-hand facet as on the 
right-hand side. We will analyze each one of the cases separately. 
1st reflection on the left-hand side 

Being < ip < (fo we can have the first two reflections on the left-hand facet, it being 

in this case guaranteed, as we assume above, that 3 or more reflections will exist. If on the 

other hand, the second reflection is on the right side, an initial part of the trajectory can 

always be represented by the first three illustrations of figure |9] (assuming < if < (po), which 

guarantee, also in this case, the existence of a third reflection B3. In order to prove what we 

have just finished saying, it will be enough to prove the ascendant nature of the sub-trajectory 
> 

-61-62- 

We establish on the parabola of the left side (illustration (a) of figure [9]) the first point of 
reflection Bi. For whatever Bi is it is always possible for us to trace an initial sub-trajectory 
BqBi with its origin in an entry angle p > ipQ. As in appendix |a] (page 14) we showed the 



sub-trajectory B1B2 which would follow it to be ascendant, the same will necessarily come 

> 

about for whatever < p < ipo may be, given that in this case BqBi will represent a more 
accentuated negative slope. Since in appendix |a (page 16) we characterized -B2-S3 only with 



basis in the ascendant nature of the sub-trajectory preceding B1B2, the conclusions to which 
> 

we arrive for -62-63 are equally valid for < 99 < (^o- 
1st reflection on the right side 

Also in this case we can have the first two reflections on the right-hand facet, it being 
guaranteed that 3 or more reflections will exist. If this does not occur, we will necessarily 
have a trajectory with the aspect of the trajectory ^0-61-62-63 illustrated in the scheme of 
figure 12, where as well there are represented two auxiliary trajectories (the dotted lines), 
^0-61^2 and ^1-62^3, which, on passing through the foci of the parabolas, present the sub- 
trajectory posterior to the reflection horizontal. IIaving_a^ a basis the laws of reflection, 
we can succinctly deduce the following: as the angle ^0-61^2 must be interior to the angle 
-60-61-62) we conclude that B1B2 is of a ascendant nature; as B1B2B3 is necessarily an interior 
angle to yli-62^3, we conclude that B3 must be situated between Ai and ^3, which guarantees 
the existence of a third reflection. Thus is the proof of the theorem [2] concluded. 
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Figure 12: Illustrative scheme to the study of the trajectory of particles with entry angle 
< (p <(po- 
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